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Abstract 

We study the competition between Josephson and charging energies in two- 



dimensional arrays of ultrasmall Josephson junctions, when the mutual ca- 
pacitance is dominant over the self-capacitance. Our calculations involve a 



VO \ combination of an analytic WKB renormalization group approach plus non- 

perturbative Quantum Monte Carlo simulations. We consider the zero frus- 

(^ ■ tration case in detail and we are able to make a successful comparison between 

^r , our results and those obtained experimentally. 



I. INTRODUCTION 
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Recent advances in submicrometer technology have made it possible to fabricate rela- 
tively large arrays of ultrasmall SIS (superconductor to insulator to superconductor) Joseph- 
son junctions [f]-f|. These arrays consist of Josephson junctions made of superconducting 
islands separated by an insulating barrier Jl|. The areas of these junctions can vary from 
a few microns to submicron sizes. In the latter case the effective capacitance of the junc- 
tions can be smaller than femtoFarads (fF=10 x 10~ 15 Farads). Under these circumstances 
the long range phase coherent properties of the Josephson junction arrays (JJA) depend 
crucially on the interplay between the Josephson (Ej) and charging energies (Eq). In this 
paper we will be concerned with studying this interplay for the specific parameter ranges of 
the experiments carried out in Delft Q. The charging energy associated with adding a single 
charge to an island is Ec e = -§q-, whereas the corresponding energy necessary to transfer a 

2 

charge from an island to a nearby one is Ec m = ^§j-- Here e is the electronic charge. For 
the Delft samples, the self capacitances are typically C s ~ 3 x 10~ 18 F, while the mutual 
capacitances are on the order of C m ~ 1 x 10~ 15 F. This means that C m is three orders of 
magnitude larger than C s . Prior to the advent of these types of JJA, most theoretical studies 



*To appear in Proceedings of the NATO Advanced Research Workshop on Mesoscopic Supercon- 
ductivity, F. Hekking and G. Schon, eds. Physica B, North Holland (1994). 

1 



assumed that the dominant term was the self-capacitive contribution, since the calculations 
are easier and also no significant changes in the results were expected from having the extra 
C m contribution. In the Delft experiments a phase diagram was obtained as a function of 
temperature vs the quantum parameter 

--£ (1) 

The phase diagram that includes the experimental results plus our theoretical results is 
shown in Fig. 1. When a = 0, the array is modeled by the two-dimensional XY model and 
it has been found to be a faithful representation |7j] of the Berezinskii-Kosterlitz-Thouless 
(BKT) scenario |§||. As a increases (in the Delft experiments the values of a cover the 
range 0.13 to 4.55, while in the Harvard samples the a can be as large as 33), the critical 
temperature decreases monotonically. As a increases further, at fixed temperature, one is 
expected to change from a superconducting phase characterized by quasi-long range phase 
coherence (SC) to a charge dominated insulating phase (I). For small and fixed a, as tem- 
perature increases, we move from a SC phase to a resistive or conducting one. For larger 
values of a the transition is expected to be between an insulating and a conducting phase. 
In the insulator normal phase boundary, the idea is that in the case when the C s = 0, a 
mutual capacitance model would map to a two-dimensional Coulomb gas model and thus 
a charge unbinding BKT transition would ensue [jT0|-|l2ll- This scenario has not been found 



experimentally since in practice, although C s can be three orders of magnitude smaller than 
C m , the electrostatic screening length is much smaller than the lattice size and the possible 
BKT transition is masked. This means that experimentally the boundary I to N is strictly 
speaking not a thermodynamic phase boundary but a cross over region |j],|J. The situation 
is different in the SC to N boundary since that is a true thermodynamic phase boundary. An 
important question is, however, if the models that have been proposed to describe the JJA 
are faithful representations of the experimental systems studied. This question is non trivial 
since there is evidence for having extraneous elements in the arrays, e.g. stray charges, that 
have not been accounted for properly in the models so far considered. It is for this reason 
that we have thought it important to carry out detailed quantitative calculations of the SC 
to N phase boundary and to make detailed comparisons with experiments. We can already 
note that our results shown in Fig.l, coming from a WKB renormalization group calculation, 
fit the experimental results quite well in the small a range. On the other hand, a rather 
good fit to the experimental results comes from carrying out a nonperturbative quantum 
Monte Carlo calculation. In this paper we describe mainly the WKB-RG calculation, and 
we will briefly mention the QMC results, which will be discussed in more detail elsewhere 
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An important result emerging from the WKB-RG calculation is the possible existence of 
a low temperature instability. This QUIT (QUantum fluctuation Induced Transition) was 
originally found in a self-capacitive model |15|JT6|| . Here we show that the results are also 
true in the case when C m dominates. This extension is non trivial since the self-capacitive 
model has no I to N phase boundary with a BKT type charge unbinding transition, even 
theoretically. One can then wonder if the competition between the insulating and supercon- 
ducting phases at sufficiently low temperatures may quench the QUIT. Within the context 
of the WKB-RG calculation this is, in fact, not the case as we discuss in this paper. From 



the WKB-RG analysis the QUIT instability may be conjectured to be a phase boundary 
from a SC to a normal phase ||15| . 

There was early support for the existence of the QUIT in the C m = 0, C s 7^ case from 
nonperturbative quantum Monte Carlo simulations |16| . In the / = case it was found that 



the helicity modulus had a discontinuity between two states with finite superfluid densities 
each. This lead to the surmise that the QUIT is in fact a transition between two superfluid 
states, one dominated by thermal and the other by quantum fluctuations. Later studies 
considered the case of / = 1/2, the fully frustrated limit, which gave a larger discontinuity in 
the superfluid density and at higher Tq UIT (f = 1/2) temperature jOJ. As we shall see from 
the WKB-RG calculation described below, the Tq UIT (f = 0) ~ Tq UIT (f = 1/2). However, 
what was found in the QMC calculations was that TquitU = 1/2) ~ 10 x TquitW = 0)- m 
a recent calculation which uses an improved QMC algorithm, that reduces the slow approach 
to equilibrium [pjj , the large hysteresis loop seen before in the / = 1/2 case was suggested 



to be directly related to the long lived metastable states seen in the classical fully frustrated 
XY model. The authors, however, did not carry out the calculations to sufficiently low 
temperature to test the / = 0, or / = 1/2 QUIT. The question then remains as to the 
physical nature of the QUIT, if it does exist. From the WKB-RG analysis there is clearly 
a low temperature QUIT instability. The known QMC results for the / = 0, C m = case 
suggest a QUIT between two different superconducting states. Further studies are needed 
to ascertain the nature of the putative QUIT. As we also mention later in this paper there 
is an experimental low temperature instability that could very well be a manifestation of 
the QUIT. 

Here we shall report on our main WKB-RG and QMC results and sketch how they were 



obtained while a fuller more detailed description will appear elsewhere JT3].[T4 . 

The outline of the rest of this paper is the following: In Sec. II we define the model 
studied here and discuss the path integral representation of the partition function used 
in our analysis. In Section III we outline our WKB-RG analysis in the case where the 
Josephson energy dominates. We obtain generalized RG equations from which we obtain 
the phase diagram, which is then compared with the experimental results 0], as well as 
with the Monte Carlo simulation results fPHf . As we mentioned above, if £e is large in the 
insulating phase the charges interact with a logarithmic potential and a BKT type transition 



must follow |T0|-[T2|. In section IV we address the question as to what happens when we 
perturb the Coulomb gas model (C s = 0) with the Josephson term in the hamiltonian; Is 
there a QUIT in the charge dominated phase? We have carried out a WKB-RG perturbative 
analysis in the insulating region and found that there is no analog to the QUIT transition 
found in the SC phase, at least to the lowest order in the perturbative analysis. 

II. THE MODEL 

A Josephson junction array can be modelled by a periodic lattice of superconducting 
islands separated by insulating barriers. Each island is characterized by a Ginzburg-Landau 
order parameter \l/(f^) = l^oO^le > where r*j is a two-dimensional vector denoting the 
position of each island. Each one of the islands becomes independently superconducting 
about the bulk transition temperature T c0 . When the temperature is lowered further, the 
magnitude of the order parameter, |\l/o(^)|) is nonfluctuating and the onset of long range 



phase coherence is responsible for the zero resistance in the arrays. The onset of phase coher- 
ence is due to the tunneling of Josephson currents between the islands and it is characterized 
by a sharp drop to the zero resistance state. The onset temperature in the arrays can be 
significantly modified by making the junctions ultrasmall. When the junction's capacitance 
is small the charging energy, i.e. the energy necessary to transfer a Cooper pair between 
the islands, can be large to the point where no Cooper pairs can tunnel any more, and thus 
the Josephson current can be reduced to the point where it is completely quenched. The 
competition between the Josephson tunneling and the charging energy, in zero field, can be 
modeled by the hamiltonian, 

H = H c + H J = ^-J2 fl (ri)C-\r 1 ,r 2 )n(r 2 ) + E J £ [1 - cos^fa) - 0(f 2 ))] (2) 

ri,r*2 <Fi,f2> 

where q = 2e. Here <j){fi) is the quantum phase operator while fi(fi) is its canonically 
conjugate number operator, which measures the excess number of Cooper pairs in the island 
placed at r*j. These canonically conjugate operators satisfy the commutation relations 

[n(n)A(r 2 )} = -iSftfr, (3) 

which imply that in the phase representation we can write ri(fj) = ~fl gjy.\ - A Josephson 
junction is essentially a capacitor with an insulating barrier. The geometric capacitance 
matrix C(?i,r 2 ) is assumed to include only the island capacitance with respect to the 
ground plane, C s , and the mutual capacitance between nearest neighbor islands, C m . In the 
square array case, of interest here, we can write the capacitance matrix as 

Cft, r 2 ) = (C. + 4C m )^, r - 2 + C m £ 6 ?if2+s , (4) 

d 

where the sum in the second term is over the nearest neighbor vectors d. The electro- 
static potential associated with this capacitance matrix model leads to the screening length 
£,e = vt?'- This means that when C m = the charges are maximally screened. In the 



arrays fabricated so far the screening length, as calculated from Eq. (^) is between 15 and 
20 plaquettes [fUJj. This means that screening effects in these arrays are important and 
that including C m in the analysis is essential. As mentioned in the introduction, however, 
prior to the fabrication of the arrays with ultrasmall capacitances, many theoretical studies 
concentrated either on mean field approaches [22], 23], which neglect the special nature of 



the BKT ordering in two- dimensions, or in 2-D studies with C m = [15,16]. Initial MC 
studies of the off diagonal problem have also appeared |1| . Since the samples fabricated 
have a typical ratio C m ~ C s 10 3 , it is clear that if one wants to make a direct quantitative 
comparison to experimental results we need to extend the analysis to include C m explicitly. 
This is the purpose of this paper. 

We should stress that the model defined by Eqs.(2-4) is not the most general model one 
could consider, for it assumes that all the electric fields are confined to the two-dimensional 
plane of the array. This may not be quite the case for real samples but one would expect that 
in the regime where £# is sufficiently large this model should be appropriate. As we shall 
discuss below the fact that one is able to make successful comparisons between experiment 



and theory in the superconducting to normal regime leads us to believe that in that regime 
at least the model gives a correct representation of the experimental system. 

There are other possibly important elements missing in the model studied in this paper. 
To wit: (i) We do not consider here the influence of a constant external magnetic field. The 
physics in that case, even in the classical limit, is highly nontrivial and requires a special 
treatment p0 |. (ii) We do not include the self- induced magnetic fields, that can also lead 
to interesting new physics [[TJJ. The reason is that the critical currents in the arrays are 
found to be too low to lead to any significant magnetic field screening effects, (iii) We do 
not include quasi-particle dissipation. Here again it has been found experimentally that 
in the SC to N phase these effects are not important ||11|. A low temperature study that 



includes quasi-particle dissipation for the diagonal capacitance matrix model has shown that 



dissipation counteracts the charging effects so as to strengthen the Josephson tunneling pi 



(iv) We could also consider the effect of disorder, in particular the one due to random stray 
charges in the array. These random frustration charges may also be more important in the 
insulating phase than in the superconducting one, which is the one we concentrate on in 
this paper. 

Here we are interested in calculating the thermodynamic properties of the model defined 
by IK. The corresponding partition function is defined by 



Z 



Tr{( 



-m\ 



(5) 



The trace is taken over <p or h. It is convenient, for calculational purposes, to change from 
the operator expression to the imaginary-time Feynmann path integral representation of Z. 
We accomplish this by using 



< n(fi)|0(f 2 ) >= <Vi,f 



exp{m(fi)(/>(fi)} 

V2^ 



(6) 



For this problem we follow steps parallel as those described in Ref. [|16j. This means that 
we discretize the imaginary-time axis into L T time slices separated by a distance e = |-. We 
recover, in principle, the fully quantum results in the limits L T —>■ oo, e — > with f3h kept 
fixed. Following Ref. |16| we get the partition function expression (up to order 0(1/L T )) 



L T -l 



fir, 



L T 



Z=U V^tjCinjf \/^#(^>) E exp 



T=0 

where we defined the action 
1 
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m(r,r)=— oo 



j:S[{4>h{m}] 



(7) 



-5[{0},{m}] 
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L T r=0 



L T 



Hj{{<f>(?,r)})-—^J2 [0(ri,rH-l)-0(f 1 ,r)+27rm(f 1 ,r)] x 



2(3q 



<r\,T2> 



xC(ri,r 2 )[0(r 2 ,r+l) -(fr(r 2 ,T) + 27rm(r 2 ,r)] 
+0(1/L T ). 



(8) 



Here Hj({<ft(r,T)}) is the Josephson hamiltonian in terms of the phase variable <f)(r,r). 
The quantum nature of this c-number functional integral comes from imposing the periodic 
boundary condition 



0(f,L T ) = 0(f,O). (9) 

These are the basic equations used in the WKB-RG and Monte Carlo calculations described 
in the next section. 



III. WKB-RG STUDY ABOUT THE SC TO N PHASE BOUNDARY 

Our approach here is to perturb the physics described by Hj({(f)(r, t)}), in the limit 
when the charging energy is small. When a = 0, Hj describes the physics of the classical 
2-D XY model |§[|. In this case we have the BKT scenario that depends on the thermal 
nucleation of vortex-antivortex pairs (VAP). The density of VAP increases exponentially as 
the temperature rises until they unbind at the critical temperature Tbkt(« = 0) = Tg^ T = 
(7rEj/2k-B)- The BKT scenario is best understood in terms of a renormalization group (RG) 
analysis @||. The RG flow diagram is obtained from a perturbation expansion in powers 

of the vortex pair density y = y§e~~ K . Here K = /3Ej and y^ is the initial condition for 
the bare vortex pair density. In the standard BKT picture there is a line of fixed points 
for < T < Tg^ T) with algebraically decaying correlation functions. In the self-capacitive 
model, at T = 0, one can map the problem to an anisotropic three-dimensional XY model, 
which must have a standard phase transition at a critical value of a c s . Around the T = 
critical point a c s) one expects to have exponentially decaying correlation functions while 
at a c s the correlations must decay algebraically. The important question is then: How do 
we go from the BKT regime, with algebraic decaying correlation functions to the very low 
temperature exponentially decaying one? There must be a discontinuity in going from one 
limit (a s = 0, T ^ 0) to the other one (a s ^ 0, T = 0). Below we discuss the evidence we 
have found, including a possible experimental candidate, that there may indeed be a QUIT 
at low temperatures. 

The situation when C s = and C m ^ is actually quite different. In that case the T = 
limit can be approximately represented by two coupled three-dimensional XY models, one 
describing the phase degrees of freedom and the other the charges. As a function of a m we 
can go from a phase dominated region, with a 3-D type XY model critical properties to one 
dominated by a 3-D Coulomb gas. There is not much known about the critical properties, 
and in particular the correlation function behavior, of the two coupled 3D XY models. So, 
strictly speaking, we can not state what kind of crossover we should expect when going 
from T = T BKT (a m ,a s ) to T = 0. Some understanding of the physics in this limit can be 
obtained by using the Villain transformation, both for the charging energy term and for the 



phase contribution | 12| . The Villain approximation is, however, valid only in a restricted 
range of a values which do not cover the full experimental range. One could conjecture, 
however, that the properties of the Villain approximated models is in the same universality 
class as the full coupled XY models, which is in fact the case when a = |]], but this needs 
to be explicitly shown. Furthermore, the general case treated here where both C s ^ and 
C m t^ is more complicated since the effective Coulomb gas in the insulating phase has a 
finite screening length. All these issues need to be studied further. 

To find the corrections to the BKT scenario due to the charging effects we carry out a 
semiclassical or WKB analysis of the model. This was originally done for the self-capacitive 



model in Ref. |15[ . Here we follow a similar approach, except that technically the problem 
is more demanding. 

To evaluate the partition function in the SC to N, or small a regime, we notice that we 
can extend the range of integration of the phases in Eq. (0) from [0,2ir] to [—00, 00], while 
at the same time all but one of the summations over the set of {m}'s can be eliminated. 
The resulting expression for the partition function, in the L T — > 00 limit, then reads 



(10) 



m(r)=— od t=0 r 

r 1 rim 
x exp — — / 

L n Jo 



where the euclidean lagrangian is 

2 



L * = \{-) £ ^rMn^Ar^ + HjiW). (11) 

\"/ <ri,r*2> 

The boundary condition, Eq. (|9|), now reads 

0(f, /3^) = 0(f, 0) + 2vrm(f), (12) 

where the {m(f)}'s are the winding numbers. We note that since the lagrangian is invariant 
under the transformation 0(f, 0) — > 0(f, 0) + 27r/(f) for all integers {^(f)}, we can extend 
the limits of integration over 0(f, 0) to [— 00, 00], the difference coming only from an overall 
multiplicative constant. Now that the limits of Eq. (|10D are all from [—00, 00] we can make 



the following change of variable [E4 



2tt — 

</>{r, f) = -7^ m (r) t + 0(f) + (p f (f, r). (13) 

Here (pf(r,r) represents the quantum fluctuations of the path about its mean value 0(f). 
These quantum fluctuations become larger than the thermal ones as a increases or as the 
temperature decreases. This means that we would need to take higher order harmonics in 
the Fourier series into account when the quantum effects are not relatively small. Because 
of the periodicity in Eq ([12]) , 0/(f, r) can be expanded in the Fourier series 

00 
Mr,r) = ((3h)- 1/2 ElM?)e^ + C.C.}, (14) 

fe=i 

where the uik = 2Trk/j3h are the Bose Matsubara frequencies. Substituting Eqs. ([13]) and 
(|i~4]) in Eq. (|T0|), expanding the Josephson term up to second order in 0fc(f), i.e. up to 
order 0(q 2 ) or equivalently to 0(a), we obtain an effective action for the classical variables 
0(f) after evaluating the integrals. In obtaining the effective action we note that once the 
integrations over the (f>k(r)'s are carried out, the partition function still includes a summation 
over the m(f)'s. In the semiclassical limit the contributions to the partition function from 
configurations with m(f) different from zero are exponentially small, so that we can safely 
take m(r) = for all f A most important property of the Josephson hamiltonian Hj is that 



it is a periodic function of its argument. This implies that in the expansion the second order 
derivative with respect to the argument in Hj is proportional to Hj itself. Specifically, for 
the cosinusoidal form of Hj we have H'j = —Hj + constant. This important property of 
Hj allows us to write the effective partition function as a 2-D classical XY model with an 
effective coupling constant. The effective partition function to this order of approximation 
isQ 

Z eff = / 1] ^ exp[-fl*l*>@)], (15) 

r 

where the effective inverse temperature is given by 

Pes = P - ^p [C-^O) - C~\d)\ . (16) 

Note that we have explicitly used the fact that C(fi, r 2 ) = C(|ri — r 2 |), valid for a periodic 
lattice, so that we can Fourier transform the capacitance matrix. 

Once we have a hamiltonian which is just like the 2-D classical XY model, we can write 
down the corresponding effective RG recursion relations to lowest order in x as 

dK , „ „ (1 - xKf , s 

H'^^W^TY (17) 

f = [2 - TK(1 - xK)]y. (18) 

In writing these equations we defined the variables 



Q 

x 



2 



VIE, 



C-\0)-C-\d) , (19) 



vr 2 7T 2 



y = eaq>[-—K(l - xK)\ = exp[-—K eff ]. (20) 

The variable x is the a parameter when the capacitance matrix is not just the self or the 
mutual capacitance. The RG equations are solved using as initial conditions K e ff(l — 0) = 
K®jj and y(l — 0) = y Q . As written, the RG equations are valid for an arbitrary ratio 
between the self and the mutual capacitances. We first notice that in the x = limit the 
RG equations reduce to the standard Kosterlitz RG equations |||9|], as they should. The 
form of the vortex density y is most important. As mentioned above, y(x = 0) grows 
exponentially with temperature. When i^O and as a function of temperature, y exhibits a 
low temperature minimum. This is shown in the discontinuous line in Fig. 2. This behavior 



1 We must note that in obtaining the effective partition function we assumed that the phases 
took values between [— oo, oo], whereas in the classical XY model the phases are constrained to lie 
in the [0, 2ir] range. Following this route makes the derivation of the effective action more direct. 
However, in the small a regime of interest here the differences between the two ranges for the 
phases can be shown to be exponentially small. 



for y is easy to understand physically. At high temperatures the difference between y and 
y is very small. However, below the minimum the increase in the vortex pair density is due 
to nucleation of VAP via quantum phase slips. Of course, we need to remember that we 
have done a perturbative calculation in x and therefore we may not be on safe ground when 
y starts increasing again at low temperatures. Nevertheless, as often happens with WKB 
derived results, the fact that the perturbative analysis shows a low temperature instability 
is likely to be true. In fact we also have found numerical evidence for the low temperature 



instability in our QMC calculations [16 



The RG equations have two nontrivial fixed points, one that corresponds to the effective 
BKT thermal fluctuations driven transition, and the other that corresponds to the QUantum 
fluctuations Induced transition (QUIT) [|15|.|TC 



The RG level curves in the (y, K) phase space, to lowest order in x, result from solving 
Eq. (|0D and 

ttxK - vrlnK + 2ir 3 y 2 = A, (21) 

K 

for different initial conditions. Figure 2 shows the RG flows for different initial conditions 
starting with different values for (yo, K°) along the discontinuous line in the figure. Each 
RG flow line corresponds to a different temperature with the arrows indicating the direction 
of increasing I. We clearly see from the figure that we can divide the temperature axis into 
three different regions. In the region between [KqIj^^Kq^], as the value of I increases 
we eliminate VAP, with the unusual property that the vortex density can initially grow 
for a while before tumbling to the critical line y = 0. This means that in this region at 
I = oo there are no VAP with infinite separation, i.e. unbounded. Below Kq^ it , as I 
increases, the RG trajectories grow away from the y = line, nonmonotonically, indicating 
that the perturbation expansion in y is no longer valid. If we associate the instability in the 
perturbation theory with the normal state behavior, one could say then that this behavior is 
characteristic of a reentrant phase transition, i.e. going from N to SC to N. This is certainly 
the case in the high temperature regime but not necessarily so at low temperatures. The 
single line that divides the two types of behaviors mentioned above is the separatrix that 
determines the critical temperature. This is the line with the highest temperature for which 
we can touch the y = line. The corresponding separatrix value of the constant A = A c is 
determined from the condition that it passes through the point {Kg KT ,y = 0). This leads 
to the result that the critical point is obtained from solving the equation, 

2 
wxK c - nlnK c - — 2vr 3 exp{-vr 2 ir c (l - xK c )} = A c {x), (22) 

K c 

where A c (x) is obtained from A c (x) = 2irxK c — it — 7rlnK c , and K c is the solution to the 
equation 2 = nxK c (l — xK c ). We'll come back to the problem of determining K^ KT in 
the next section, where we make a comparison with the experimental results. It is known, 
however, that the determination of Tbkt using the RG equations is not quantitatively 
exact, since the RG analysis is explicitly derived for the Villain action. In comparing with 
experiment in the next section we will take this into account. Here we present the corrections 
to the classical results to the leading order in x, which give the correct qualitative trends. 
Specifically, expanding in powers of x we find that Tbkt and Tqjjjt are given by 

9 



T BKT ^T^ T -^x + 0(x 2 ), (23) 

fc B 

T QmT ^^x + 0(x 2 ). (24) 

Kb 

Notice that these equations are applicable not only in 2-D, for if the system described 
by Eq. (|15|) has a transition point at some K^ s then the equation K^ s = K — xK 2 has two 
solutions for K, which are the ones implied in Eq. fll8l). Moreover, notice that the results 
to the first order in x are independent of the specific value of T B ^ T . This means that if 
we consider the finite magnetic field case, the corresponding critical temperature will be 
T C (B) « T(°)(B) — {Ej/k-Q)x + 0(x 2 ). Furthermore, we notice that, to the lowest order in 
x, the Tquit must be the same with and without a field. This fact will be compared with 
the experiments in the next section. 

The explicit leading order calculation of the correction to the BKT critical temperature 
in the asymptotic limits in which either the self or the mutual capacitance dominates results 
in 



JBKT _ I 3tt Ej ^ W 

T (0) I 1 4 E Cm n 



(Eg, 

\Ej 



( E Cm 
V Ej 



if c s > c m 

(25) 
if C s < C m 



here z is the coordination number of the array, and for a square array in two dimensions, 

z = A. 



IV. COMPARISON TO EXPERIMENT 

We now move on to a brief discussion of how we obtained the results presented in Fig. 
1. As mentioned before, in trying to find quantitative correspondence between experiment 
and theory it is important to ascertain the validity of the theoretical models employed to 
study the arrays. We have carried out two different types of checks. One based on the RG 
analysis described in the previous sections and the other from a nonperturbative quantum 
Monte Carlo calculation |l3 , [l4l . We discuss the RG analysis here and only briefly mention 



the QMC results, with more details left for forthcoming publications [0,0 • 

As mentioned in section III , it is known that the RG equations do not lead to quanti- 
tatively exact results for the critical temperatures. They do, of course, lead to the correct 
universal critical exponents. However, what has so far been measured experimentally is 
the phase boundary between the SC and the N phases. We then need a consistent way to 
compare the RG results with the experimental results. 

We first note that the phase diagram of Fig. 1 is plotted as a function of a m = -f 2 ^, since 
experimentally a s is three orders of magnitude smaller. We can then write K e ff = K(l—xK) 
with Cs = as 

K eff = K- ^K 2 . (26) 

Next we set the critical temperature for the classical model to be the one obtained in 
classical MC simulations (e.g. pOfl ) 1/A^°) w 0.93, so that the critical temperature T c (a m ) 



of the actual model is given by the equation 

10 



k B TJa r 



2id 0) 




2Om K (0) 



(27) 



To lowest order in a m this equation gives 



O^m 



T c (a) = T c (0) - ^, (28) 

whereas the maximum value of a m for which there is a physical solution is a m = |/cbT c (0) ~ 
1.4. The results obtained from this analysis are shown as a discontinuous line in Fig. 1. By 
following this approach we see that for a m < 1 the RG result is actually quite good when 
compared to the experimental and the MC results. 

We also have extended our previous QMC calculations to the case when the off-diagonal 
capacitance is dominant. The results are shown in Fig. 1 by the crosses, including their error 
bars. It is clear from these results that the correspondence between experiment and QMC 
results is excellent, up to nonperturbative values of a m . This leads us to the conclusion that 
the model studied here does provide a good representation of the experimental system, at 
least in the SC to N regime. In the next section we will briefly discuss what happens in the 
insulating region to normal region. 

The discussion presented above dealt with the SC to N phase boundary. What about 
evidence for a QUIT? In this regard we note that in the experimental results of Ref. || there 
are results for a sample with a nominal a m = 1.67 for which there is a double type of reentrant 
behavior. The low temperature glitch seen in the resistance versus temperature diagram 
occurs at T" = AOmK. Moreover the latter instability is also characterized by an increase in 
noise fluctuations in the IV characteristics measured in this sample. If we assume that what 
is seen at T" is related to the QUIT, using Eq.(24) and the parameter values of the experiment 
we get a Tquit = 33mK, rather close to the experimental value. Furthermore when the 
same experiment is repeated in a small magnetic field the low temperature instability is 
found at he same temperature, i.e.; T'(f = 0) = T'(f = 0.08). This result is also consistent 
with Eq. (24) which also leads to a Tquit independent of / at leading order in a m . These 
results may just be coincidental and more work needs to be done to conclusively connect 
the Tquit with the low temperature instability already seen in the Delft experiments. 

V. INSULATING TO NORMAL CROSS OVER 

As mentioned in the introduction, as a m increases there is a SC to I transition at fi- 
nite temperature. In the a s = case, the insulating phase has been modelled as a two- 
dimensional Coulomb gas of charges with a possible BKT charge unbinding transition 
TOj-12 . This situation has been studied extensively, in particular in Ref. [12]. Here we 



ask if there is an equivalent QUIT in the insulating phase at low temperatures. If we draw 
an analogy to the quantum induced vortices in the SC phase one could also imagine that 
the the number of free dipoles in the arrays could increase due to quantum fluctuations. 
However, as we show below the vortices and charges are not dual to each other in that 
sense. 

The calculation described here aims at finding the leading correction to the charging 
hamiltonian due to Josephson junction fluctuations. We then expand the Josephson contri- 
bution to Z as 
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As before, we use Eqs. (12) and (13), but this time we integrate out both </>/(r, r) and 0(f), 
which leaves us with an effective action for the {m}'s, 



ZnZtU £ exp 

r m(r)=—oo 






m(r + d) — m(r) 



(30) 



where we have assumed that C s <C C m . The function Z^ does not contain the {m} variables, 
and the effective coupling constant K e g is 



Ki = K- 1 



2tt 2 



a,, 



g{K) 



(31) 



given as a function of K = (PE Cm )/2n. As in the vortex dominated case, we have ended 
up with a Coulomb gas problem but with a renormalized coupling constant. The function 
g{K) determines the importance of the zero point fluctuations of the phases on the charge 
dominated phase. The function g(K) is defined as 



g{K) 



1/2 



dt[l — cos(27r£)] exp 



(2tt) 



-Kt(l - t) 



(32) 



For a BKT type phase transition we have K e fj = K^ KT = 2 /n (T^j. It is important here to 
see if, within this approximation, the system shows a QUIT or a reentrant transition. We 
then need to study the number of solutions to this equation. From the fact that K^ KT ~ 0(1) 
we can see that for any value of l/a m the function g(K) has only one solution for K c . This 
means that, to this order of approximation, there is no QUIT in the charge dominated 
phase in this model. 

A question that immediately arises is: Why is there a difference between the vortex 
and charge dominated phases, in particular in view of the duality between the two phases 
extensively studied in |I2| ? The reason is that the duality is not exact since there is a 
term in the action, obtained using the Villain approximation, that breaks this symmetry. 
If one includes this term we then see that the cost of producing quantum fluctuations in 
the vortices is bounded from above whereas the corresponding cost in the charge dominated 
phase is unbounded. In our calculation we have kept these contributions intact. 
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FIGURES 

FIG. 1. Temperature vs. charging energy phase diagram. The experimental results are 
denoted by squares. The RG results are given by the discontinuous line while the quantum Monte 
Carlo results are given by the crosses and joined by a continuous line as a guide to the eye. The 
latter results include the statistical error bars in the calculations. 

FIG. 2. Renormalization group flow diagram. The discontinuous line indicates the vortex pair 
density as a function of temperature. See text for a discussion of the analysis of this diagram. 
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